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ABSTRACT 

We analyze several aspects of the recently noted neutron star collapse instability 
in close binary systems. We utilize (3+1) dimensional and spherical numerical general 
relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of 
this instability. We derive the modified conditions of hydrostatic equilibrium for the stars 
in the curved space of quasi-static orbits. We examine the sensitivity of the instability 
to the neutron star mass and equation of state. We also estimate limits to the possible 
interior heating and associated neutrino luminosity which could be generated as the stars 
gradually compress prior to collapse. We show that the radiative loss in neutrinos from 
this heating could exceed the power radiated in gravity waves for several hours prior 
to collapse. The possibility that the radiation neutrinos could produce gamma-ray (or 
other electromagnetic) burst phenomena is also discussed. 


Subject headings: neutron stars: structure — neutron stars: binaries - general relativity 
— gamma-ray bursts 
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1. Introduction 

In recent numerical studies of the relativistic hy¬ 
drodynamics of close neutron star binaries in three 
spatial dimensions (Wilson & Mathews 1995; Wilson, 
Mathews & Marronetti 1996, henceforth WMM96), 
it was noted that as the stars approach coalescence 
they appear to experience a collapse instability. For 
an appropriate equation of state, binary neutron stars 
might generally become black holes many seconds 
prior to merger. If correct, this effect could have a sig¬ 
nificant impact on the anticipated gravity wave signal 
from neutron star binaries near coalescence. Such pre¬ 
merger collapse might also be associated with heat¬ 
ing, neutrino production, and electromagnetic bursts 
as the released gravitational energy from collapse is 
converted into thermal energy of the stars. 

Moreover, the numerical evidence that such an in¬ 
stability exists poses a number of new questions such 
as the sensitivity of the instability to the specific 
equation of state employed, or the intrinsic spin and 
masses of the stars. One would also like to understand 
the time history of the collapse and any associated 
electromagnetic or neutrino emission. 

In this paper we present some new three dimen¬ 
sional (3D) calculations which begin to examine these 
issues. Unfortunately, however, such relativistic hy¬ 
drodynamic calculations in three spatial dimensions 
are computationally expensive. A complete system¬ 
atic study of this instability in three spatial dimen¬ 
sions will be long in coming. In this paper, however, 
we show that in large part this effect can be replicated 
in terms of modified one-dimensional spherical rela¬ 
tivistic hydrodynamics. We show that the relativistic 
effects of placing stars in a close binary can be ap¬ 
proximated by adding a term involving an average 
Lorentz-like factor which increases the effective grav¬ 
itational forces on the stars. The collapse observed 
in the three dimensional calculations can be under¬ 
stood in this one-dimensional framework and one can 
survey easily the sensitivity of this effect to param¬ 
eters characterizing the binary and the neutron star 
equation of state. 

We can also follow the possible precollapse com¬ 
pression and heating of the neutron star material. 
This provides a framework in which to model the 
possible associated neutrino and electromagnetic sig¬ 
nals such as gamma-ray bursts. We show that sig¬ 
nificant heating and neutrino emission is possible as 
the stars gradually compress before they reach the 


collapse instability. During the heating epoch the as¬ 
sociated neutrino and electromagnetic radiative losses 
may dominate over the power loss from gravitational 
radiation. 

2. Field Equations 

Our method of solving the field equations in three 
spatial dimensions was discussed in Wilson & Math¬ 
ews (1995) and WMM96. Here we present a brief 
review of some features relevant to the present dis¬ 
cussion. 

We start with the slicing of spacetime into a one- 
parameter family of hypersurfaces separated by dif¬ 
ferential displacements in time-like coordinates as de¬ 
fined in the (3-1-1) formalism (Arnowitt, Deser & Mis- 
ner 1962; York 1979). 

Utilizing Cartesian x, y, z isotropic coordinates, 
proper distance is expressed 

ds^ = —(of^ — I3il3^)dt^ + 2f3idx^dt + (j)‘^Sijdx^dx^ (1) 

where the lapse function a describes the differential 
lapse of proper time between two hypersurfaces. The 
quantity fii is the shift vector denoting the shift in 
space-like coordinates between hypersurfaces. For an 
orbiting binary, is dominated by the orbital motion 
of the system plus a small contribution from frame 
drag (WMM96). In the frame of one star in the one 
dimensional calculations described here, most of the 
effect of the (di is transformed away. 

The curvature of the 3-geometry is described by a 
position dependent conformal factor <j>^ times a flat- 
space Kronecker delta. We refer to this gauge choice 
as the Conformally Flat Condition or CFG. For a 
static system, the vanishing of the Weyl tensor in 
three dimensions guarantees (cf. Weinberg 1972) that 
there exists a conformally flat solution to the Ein¬ 
stein equations. One must be careful, however, not 
to overimpose symmetry conditions (e.g. Cook et al. 
1996). For a dynamic system, one can always im¬ 
pose conformal flatness as an initial condition. There 
are, however, nonzero time derivatives of the spatial 
metric and extrinsic curvature which can begin to in¬ 
troduce off-diagonal elements of 7 ^ as the system 
evolves. In particular, the imposition of conformal 
flatness excludes gravity wave information contained 
in the transverse traceless components of the metric. 
However, as discussed below, several recent studies 
have indicated that this approach is still an excellent 
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approximation when a comparison with exact results 
can be made. 

The implementation of the CFG means that we 
solve the constraint equations of general relativity at 
each time as though there were a fixed distribution 
of matter. We then evolve the hydrodynamic equa¬ 
tions to the next time step against this background 
metric. Thus, at each time slice we can obtain a solu¬ 
tion to the relativistic field equations and information 
on the hydrodynamic evolution. Information on the 
generation of gravitational radiation can then be ob¬ 
tained from a multipole expansion (Thorne 1980) of 
the transverse traceless components of the metric. 

It is important to appreciate that at each time slice 
a numerically valid solution to the field equations is 
obtained. In this way the strong field properties of the 
system are included. For this reason, this approach is 
a significant improvement over a post Newtonian ap¬ 
proach (which is also conformally flat at low order: cf. 
Appendix A). The hydrodynamic variables respond to 
these fields. An approximation we have made herein 
is the neglect of an explicit coupling of the gravity 
waves. These, however, contribute negligibly to the 
metric, stress energy tensor, or hydrodynamic evolu¬ 
tion (WMM96). When desired, this coupling can be 
added via a multipole expansion. 

We reduce the solution of the equations for the field 
variables (j>, and a to simple Poisson-like equations in 
flat space. We begin with the Hamiltonian constraint 
equation (York 1979) which reduces to (Evans 1985; 
WMM96), 

V^<f> =-47rpi. ( 2 ) 

In the Newtonian limit the source term is dominated 
by the proper matter density p. In the strong field 
of the orbiting binary, however, p must be enhanced 
by a generalized curved-space Lorentz factor W [cf. 
Eq. (11)]. This derives directly from the occurrence 
of the four velocity in the stress energy tensor. There 
are also contributions from the internal energy density 
ep, pressure P, and extrinsic curvature Kij. Thus we 
write. 


Pi = TT 


pW^ + pe(^rW^ - T + 


- 

IbTT 


(3) 


where F is an adiabatic index from the equation of 
state as defined below. Similarly, the lapse function 


is determined from, 

(«()))= 47r/92 , (4) 


P2 = 


a(j)^ 


p{3W^ -2)+ pe[3F(IF2 -b 1) - 5] 




(5) 


In WMM96 it was shown that the collapse insta¬ 
bility can at least in part be traced to the effect of the 
Lorentz-like factor W on the source density for <j) and 
In WMM96 and below it is shown that terms 
which scale as (W^ — 1) also enter into the hydrody¬ 
namic equations in a way which enhances the gravita¬ 
tional force on each star. In Appendix A we suggest 
that, even in a post-Newtonian approximation, such 
terms might cause the effective gravitational poten¬ 
tial to be deeper than it would be for static isolated 
stars. 

Regarding the reliability of the CFC as an ap¬ 
proach to this problem we note that a recent study 
(Cook, Shapiro & Teukolski 1996) has shown than an 
axially symmetric CFC approximation is quite good 
when computed physical observables are compared 
with the exact results for axisymmetric extremely 
rapidly rotating neutron stars. This is the simplest 
system for which an exact metric begins to differ from 
a CFC metric. 

In another recent application to the nonaxisym- 
metric case of orbiting binaries Reith & Schafer (1996) 
have shown that an expansion using this metric is 
identical to a post-Newtonian expansion for terms of 
order (v/c)'^. The first deviation appears in terms of 
order (v/c)^. However, we find that the deviations are 
small. The expressions in their paper are in terms of 
a dimensionless parameter v = 17111712 /(mi +m 2 )‘^- It 
is common in post-Newtonian expansions to compare 
with a Schwarzschild orbit for which i/ = 0. However, 
for neutron star binaries v = 0.25 is most appropriate. 
For equal-mass neutron-star binaries v = Xj4 exactly. 
Even for mi/m 2 = 2 (a relatively large asymmetry 
for neutron stars) v = 0.22. For v = 0.25 the differ¬ 
ence between the conformally flat and post Newtonian 
{v/c)^ correction for the perihelion advance is about 
is 4.5%. The (u/c)^ term for the momentum differs 
by about 2.8%, and the angular momentum term dif¬ 
fers by 24.1%. Since (u/c)^ ~ 10“^ for the binaries 
considered here, these differences in the two-body dy¬ 
namics are probably insignificant. Note also, that in 
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the present application we compute an exact instan¬ 
taneous numerical solution to the Einstein equations 
using this metric and do not rely upon an expansion 
which may deviate in individual terms but still pro¬ 
vide the correct results. Note also, that the principle 
effect we are investigating here (that due to the W^ — 1 
terms) is of order (w/c)^ (see Appendix A) for which 
the post-Newtonian and conformally flat terms agree 
exactly. Also, note that the effect described here is 
a relativistic effect which completely dominates [see 
below] over the possible stablizing influence of New¬ 
tonian tidal distortion as proposed in Lai (1996). The 
effect we describe was not considered in that paper. 


3. Relativistic Hydrodynamics 


To solve for the fluid motions of the system in 
curved spacetime it is convenient to use an Eulerian 
fluid description (Wilson 1979). We begin with the 
perfect fluid stress-energy tensor, which in covariant 
form can be written. 


— {p + pt + P)U^U^-\-Pg^v , ( 6 ) 


where e is the internal energy per gram, P is the pres¬ 
sure, and Ui, is the four velocity. 

By introducing a set of Lorentz contracted state 
variables it is possible to write the relativistic hy¬ 
drodynamic equations in a form which is reminiscent 
of their Newtonian counterparts. The hydrodynamic 
state variables are: The coordinate baryon mass den¬ 
sity, 

D = Wp , (7) 

the coordinate internal energy density, 

E = Wpe , ( 8 ) 


the spatial three velocity, 


W = a 


Ui 

(ji^W 




and the coordinate momentum density. 


(9) 


S^ = {D + TE)U^ . ( 10 ) 


The Lorentz-like factor W is 


W = aU* 


1 -b 



UF 


2-11/2 


( 11 ) 


and the EOS index T is 

r = l + — . (12) 


Note that in flat space {a = (j) = 1),W reduces to the 
usual special-relativistic Lorentz factor. 

In terms of these state variables, the hydrodynamic 
equations are as follows: The equation for the conser¬ 
vation of baryon number {pU^).i, = 0 takes the form. 


dD _ ^^ dXogfj) 
dt ~ dt 




(13) 


The internal energy equation is derived from = 

0 , 


dE 




-p 


dW 


1 d 
(j)^ dxi 


( 0 ®ITE^') 


(14) 


The spatial components of the momentum conserva¬ 
tion condition {T^ = 0 ) takes the form. 


dS, 

dt 


-6S, 


d\og4> 


1 d 


dP 


(j)^ dxi dx^ 


2a(D + r£)(W'-T,2|S^ + 

ir(D + rE)g. 


1 dx^ 
(15) 


where for the present stability analysis we have set 
the radiation reaction term (WMM96) to zero. 

Regarding the stability of our treatment of nu¬ 
merical relativistic hydrodynamics, we have applied 
a number of standard tests (e.g. shock tubes, pres¬ 
sureless collapse, etc.) as noted in WMM96. One 
important test in the present context is that of sta¬ 
ble stars in a stable orbit (with no radiation reaction). 
We have found that for such systems, equilibrium con¬ 
figurations are obtained after a small fraction of an 
orbit. When the velocity damping is removed there 
is no discernible change of the stars for several orbit 
periods. This illustrates an advantage of the shifted 
grid which we employ. There is essentially no matter 
motion with respect to the grid once a stable equilib¬ 
rium has been achieved. Hence, numerical stability 
can be maintained for a long time. 

Another important test is that of a single non¬ 
rotating star on the three-dimensional spatial grid. 
We find that the equilibrium gravitational mass as 
a function of central density agrees with the spheri¬ 
cal hydrostatic Tolman-Oppenheimer-Volkoff equilib¬ 
rium gravitational mass as a function of central den¬ 
sity to within a fraction of a percent. We also find 
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that a dynamical instability ensues once the stellar 
mass exceeds the maximum hydrostatic mass as ex¬ 
pected. 

In isotropic coordinates, the condition of hydro¬ 
static equilibrium for the stars {dS^/dt = 0, = 0, 

d log (j)/dt = 0 ) can be inferred from equation (p^), 


dP 

dx^ 


= -(p + per) 


9 log a 
dx'- 


Uj df3^ 
a dx^ 


9 log a 
dx^ 


^ d\og(j) 

dx^ 





.(16) 


Some discussion of the relative magnitude of the 
terms in Eq. is useful. The first term with 
dloga/dx^ is the relativistic analog of the Newto¬ 
nian gravitational force. In the Newtonian limit 
a ^ 1 — GM/r. Hence —dloga/dx^ —> GM/r'^. 
In Eq. (16) there are two ways in which the effec¬ 
tive gravitational force increases as W exceeds unity. 
One is that the matter contribution to the source den¬ 
sity for a 01 (j) is increased by factors of ~ [cf. 
Eqs. (^,||)]. The more dominant effect is that from 
the terms in Eq. which scale as iW^ — 1). These 
terms result from the affine connection terms 


in the covariant differentiation of These terms 

have no Newtonian analog but describe a general rel¬ 
ativistic increase in the curvature gravitational force 
as the specific kinetic energy of the system increases. 
This increase in effective gravity as the stars approach 
each other can be thought of as a correction to the 
Newtonian gravity which scales as {W^ — 1) times 
the Newtonian gravity. This (VE^ — 1) factor can be 
thought of as a kind of specific kinetic energy [cf. Eq. 
m from the orbital motion of the binary. The extra 
dlog(j)/dx'^ term further increases the effect by a factor 
of 2. This factor comes from (jp ~ (1/a). A further 
increase of binding arises from the Kij terms in 
the field sources, but these terms are much smaller 
than the — 1 contributions. 


In our shifted spatial grid, the fluid three velocities 
E® are nearly zero. Hence we can use Eq. (^ to find 
Ui as a function of /9b We can also replace /9® in 
Eq. (^) with the dominant contribution from orbital 
motion, (d k, Rxui, where R is the coordinate distance 
of the stars from the center of mass. Hence, Eqs. (^) 
and (0) give, 


{W^) 


1 


(1 - 

1 

(l-u2/c2) > 


( 17 ) 


where c = < 1 is the coordinate light speed. 

In our simulations, an effective velocity of [loRjc) ~ 
0.25 is obtained for the last stable orbit of 1.45 Mq 
stars. In the 3D calculations the average (IT^ —1) typ¬ 
ically rises up to ^ 5% before the orbit becomes dy¬ 
namically unstable. Thus, we estimate that before or¬ 
bit instability, the effective hydrostatic gravitational 
force on the stars is increased by ^ 10 % over that of 
stationary non-orbiting stars for which (lE^ — 1 ) = 0 . 
This increased gravitational force increases the cen¬ 
tral densities as the stars approach and can induce 
collapse. 

It is of interest to compare the magnitude of the 
{W^ — 1) correction to the Newtonian gravity with the 
magnitude of the Newtonian tidal energy which would 
tend to stabilize the star. In Lai (1996) it was esti¬ 
mated that the Newtonian tidal energy should scale 
as 

AEudal « -A ^^ 6 ^" (18) 

where for neutron stars A 0.1, M is the mass of a 
star, R is the neutron star radius, and r is the orbital 
separation. In contrast, the correction to the New¬ 
tonian self gravity from the motion of the stars in a 
binary is 

GM'^ 

AEgr « 2 (W^ - 1 )^ (19) 

Taking the ratio of these two contributions, we find. 


AEudai ^ A /Ry 
AEgr ^ 2(IE2 -l)[rj 


( 20 ) 


where we have used typical values near collapse (cf. 
Table 1) of R/r ~ 0.2 and {W‘^ — 1) ~ 0.05. Thus, the 
effect of the increased relativistic gravitational force 
is expected to dominate over the stabilizing tidal dis¬ 
tortion by about 4 orders of magnitude. In a similar 
analysis, we estimate that even for white dwarfs near 
merger, this relativistic increase in the gravitational 
energy exceeds the stabilizing Newtonian tidal energy. 

The centrifugal term in Eq. (16) is dominated by 
the contribution from orbital motion Uj{dj3^/dx^) « 
UjUj ~ 10“®. This term varies little over an indi¬ 
vidual star and inside a star this term is small com¬ 
pared to the 91oga/i9a;* term. Hence, this term can 
be neglected in the discussions of stellar stability. It 
is important, however, for determining the orbits and 
gravity wave frequency (WMM96). 
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4. Equivalent Spherical Model 


The hydrodynamic equations become: 


To better understand the relativistic effects de¬ 
scribed herein, it is useful to reduce the hydrodynamic 
equations to an approximate spherical model. Such a 
model can also be used to make a schematic survey of 
the sensitivity of collapse to EOS parameters and pos¬ 
sible interior heating. From Eq. ( 0 ) we see that the 
configuration of each star can be described by a mod¬ 
ified version of the familiar equation of hydrostatic 
equilibrium. This is true as long as the contribution 
of orbital motion to (W^ — 1 ) can be treated as a con¬ 
stant factor and the KijK^^ and centrifugal terms can 
be ignored. Of course, is not constant over the 
stars. However, in our three dimensional calculations 
it is observed to vary little over the volume of a star. 

Contours of constant — from a three-dimensional 

calculation are shown in Fig. ^ From this we de¬ 
duce that it is not a bad approximation to replace 
(pj /2 _ 2 ^^ source equations and hydrodynami- 

cal equations with 

-1 + {W^-1)) , (21) 

where Wr is the contribution from radial motion in¬ 
side a star, and (lEg — 1 ) is an approximately con¬ 
stant factor which accounts for the influence of orbital 
motion in the curved space-time of the binary. The 
equilibrium and stability of a binary star can then be 
approximated using a one-dimensional description. 

However, since the metric variables a and (j) de¬ 
pend upon the density distribution, it is not possible 
to directly compute the hydrostatic equilibrium of the 
star. Instead, the star must be evolved hydrodynam- 
ically (with damping) as (Wq — 1 ) is increased to ob¬ 
tain the new equilibrium configuration. Hence, we 
construct a modified spherical hydrodynamic model 
as follows: 

For a given distribution of mass and energy, the 
Poisson equations (|^) and (^ for tp and a can be 
integrated directly. The only difference is that the 
source terms now become 

Pi « Y (p(l + e + er(IE,^ - 1 + - 1 ))^ , ( 22 ) 

and 

P2 « ^(ip(l + eT){W^-l + {W^-l)) 

-\-p{l + 6 + 6e(r — 1)) ) • 


dD _ ^^ dlogcj) 
dt ~ dt 


1 d 

(j)Gj-2 Qj. 




( 24 ) 


The equation for internal energy conservation be¬ 
comes, 




-P 


dWr 


1 d 


dt dr 

The momentum equation is 
dlogp 




(.25) 




= —6Sr 


dt 


1 d Q 2c T/r\ 

^^(0 r S,V ) - 


- {W^ + {W^-l))iD + TE)^ , (26) 


where we have neglected the centrifugal term as noted 
above. 

The calculations reported here were performed on 
an Eulerian grid in which the star is resolved into 
about 100 radial zones. 


5. Equation of State 

A key part of the calculations presented here is the 
use a realistic neutron star equation of state (EOS). 
The orbital calculations presented in WMM96 used 
the zero temperature, zero neutrino chemical poten¬ 
tial EOS from the supernova numerical model of Wil¬ 
son & Mayle (1993), Mayle, Tavani & Wilson (1993). 
Calculations made with this EOS for a model of su¬ 
pernova 1987A give an explosion energy of 1.5 x 10®^ 
ergs, consistent with observation. Also, the neu¬ 
trino spectra and time of neutrino emission are in 
good agreement with the IMB (Bionata et al. 1987) 
and Kamiokande neutrino detections (Hirata et al. 
1987). These models also reproduce the desired abun¬ 
dance distribution of r-process heavy elements in the 
baryon wind from the proto-neutron star (Woosley et 
al. 1994). The maximum mass neutron star for this 
EOS as converted to a zero temperature version for 
our present studies is Me = 1.70 M©. An important 
point is that with an EOS which would allow a higher 
mass neutron star, Wilson & Mayle (1993) were not 
able to obtain satisfactory results. 
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In WMM96 only cold equilibrium configurations 
were computed. However, in the present work we wish 
to examine the possible heating of the stars as they 
collapse. Hence, we include finite temperature effects 
in the EOS. The electron fraction is small for neutron 
stars Ye << 1. Hence, for the heating calculations 
of interest here, we can relate the temperature to the 
internal energy by assuming a nonrelativistic fermi 
gas of neutrons. 

We wish to analyze the sensitivity of the collapse 
instability to the neutron star equation of state. To do 
this we diminish the maximum mass achievable for a 
given equation of state by imposing a maximum value 
for the index T at high density. From this maximum, 
the adiabatic index asymptotes to 2 at high density 
to guarantee causality. 

We find a maximum neutron star mass of Me = 
1.55, 1.64, and 1.70 for = 2.297, 2.346, 2.470, 
respectively. This range of masses is consistent with 
(and even slightly above) the upper range of the 
observed upper mass limit for neutron stars. Finn 
(1994) has assigned a lower limit of 1.15 to 1.35 Mq 
and an upper limit of 1.44 to 1.50 Mq at the la (68%) 
confidence level. At the 2a (95%) confidence level the 
upper limit increases to 1.43 to 1.64 M©. In an in¬ 
dependent approach, Bethe and Brown (1995) have 
recently argued from nucleosynthesis constraints that 
the maximum neutron star mass is 1.56 M©. They 
also point out that if kaon condensation is taken into 
account the critical mass may only be 1.50 Mq. If 
the maximum observed stellar mass were as low as 
the Icr upper limit, i.e. 1.50 M 0 , it could be that 
almost all neutron star binaries would collapse before 
coalescence. 

With the present state of knowledge of the nu¬ 
clear equation of state at high density, however, it 
is still possible that the maximum neutron star mass 
could be significantly greater than 1.70 Mq. That 
is, the observed low mass limits may only be an arti¬ 
fact of the way in which neutron stars are formed in 
type H supernovae rather than a limit from the EOS. 
We have made some preliminary studies of stars with 
Mg = 1.45 Mq for an EOS with Me = 1.85 Mg. We 
have not observed collapse before the orbit instabil¬ 
ity is reached. It seems likely that for a sufficiently 
stiff EOS that merger will occur as two neutron stars 
as considered in many Newtonian and post Newto¬ 
nian simulations (e.g. Rasio & Shapiro 1992; Zhuge, 
Centrella, McMillan 1994; Janka & Ruffert 1996). 


6. Summary of 3D results 

6.1. Summary of Previous Results 

In WMM96 orbit calculations were made for two 
1.70 Mg baryonic mass neutron stars for an EOS 
for which the gravitational mass in isolation was 
Mg = 1.45 Mg and the critical neutron star mass 
was Mg = 1.70 Mg. Orbit solutions were sought 
for three separate values of angular momentum. The 
neutron stars were taken to be corotating initially, 
although it was noted that relativistic effects subse¬ 
quently induce some fluid motion in the stars relative 
to the corotating frame. In the present study we have 
not considered the more realistic possibilities of initial 
neutron-star spins or non equal masses. Such systems 
contain less symmetry and require a larger computa¬ 
tional effort which we leave for future work. 

The first calculation was made with an orbital an¬ 
gular momentum of 2.2 x 10^^ cm^. The stars settled 
down into what appeared at first as a stable orbit, but 
later (less than one complete orbit) the stars began to 
slowly spiral in. The next calculation was made with 
an angular momentum of 2.3 x 10^^ cm^ for which 
the orbit appeared stable. However, after about 1 to 2 
revolutions the central densities were noticed to be ris¬ 
ing. By the end of the calculation the central baryonic 
densities had continuously risen to about 2.7 x 10^® 
g cm“^ (« 10 times nuclear matter density) which is 
near the maximum density for a stable neutron star 
for an EOS with Mg = 1.70. It appears that neutron 
stars of this mass range and the adopted equation of 
state may continue to collapse as long as the released 
gravitational energy can be dissipated. For this orbit 
the stars are at a separation distance of dp/m = 9.5, 
far from merging. By the time the calculation was 
ended, the minimum a had diminished to 0.379 and 
(jp' risen to 2.05 corresponding to a minimum coordi¬ 
nate light speed of 0.18. 

A third calculation was made with the angular mo¬ 
mentum increased to 2.7 x 10^^ cm^. As can be seen in 
Tabled and the stars at this separation dp/m = 12.4 
seemed both stable and in a stable orbit. However, 
with only a slight increase in baryonic mass (Mf, 
1.598 ^ 1.620 a collapse ensues. 

6.2. New 3D Results 

In the present work, these results are supplemented 
with additional 3D calculations for initially corotat¬ 
ing stars. The new results are obtained with better 
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resolution and an improved treatment of boundary 
conditions. These calculations have been run several 
times longer than in WMM96 so that for cases where 
the instability occurs we have followed the collapse to 
higher densities and stronger fields. 

These new results along with the previous results 
are summarized in the Tables ^ through For the 
star with Mg « 1.45 Mq in isolation, we have added 
calculations at intermediate orbital angular momenta 
of 2.5 X 10^^ and 2.6 x 10^^ cm^ here we have run 
for much longer times and further into the collapse. 
We find that even at 2.6 x 10^^ cm^ the stars collapse 
while still at a distance of over 50 km apart and with 
a gravity wave frequency of only 250 Hz. The collapse 
instability appears to onset between 2.6 and 2.7 x 10^^ 
cm^ as — 0.05 (depending upon the baryon 

mass). 

The rate of energy and angular momentum loss 
generally increases as the compression proceeds and 
a becomes small. We note, however, that for un¬ 
stable stars or orbits, the systems are no longer in 
quasi-equilibrium orbits. Since the radiated energy 
and momentum are sensitive functions of separation 
and w, the computed values of energy and angular 
momentum loss become unreliable as an average es¬ 
timate. This is the reason for the lack of monotonic¬ 
ity in the E values of Tables 0-1 Eventually the 
system approaches two black holes and is no longer 
well describable in our framework. From the rates of 
energy and angular momentum loss for these calcu¬ 
lations we make a crude estimate (WMM96) that a 
delay of about 5 seconds occurs between the collapse 
instability and the orbit instability for stars with this 
EOS and Mq = 1-45 M©. This would have interesting 
consequences on the gravity wave or electromagnetic 
signals as discussed below. 

We have also run 3-dimensional calculations for 
stars which would have Mq = 1.40 Mq and 1.35 Mq 
in isolation (Tables ^ and p. For these systems the 
EOS was varied as well as the angular momentum. 
Results from these calculations are summarized in Ta¬ 
bles 2-4. We found as expected that the 1.40 Mq 
stars are stable for lower angular momenta than the 
1.45 Mq stars. However, stars of this mass will still 
collapse for an EOS which is more soft. For example, 
the 1.40 Mq star will collapse for J = 2.6x10^^ cm^ if 
the upper mass limit from the EOS is reduced to 1.55 
Mq. We have found that with the Mq = 1.70 Mq 
EOS it is necessary to reduce the mass to Mq = 1.35 
Mq before the stars can survive to final orbit plunge 


without collapsing first (see Table 3). 

6.3. Connection to ID Calcnlation 

One of the concerns in the 3D calculations of 
WMM96 was whether the stars and field variables 
were sufficiently resolved to produce reliable numeri¬ 
cal results. In the 3D calculations, the spatial resolu¬ 
tion only provided ~ 10 — 15 zones in radius across a 
star. In the ID calculations, however, one can easily 
provide many radial zones. We have made a survey 
of various physical quantities such as the central ex¬ 
trema in a, (j), and p, as well as the total gravitational 
mass Mq as a function of the number of radial zones. 
We have found that there is no significant difference in 
the field variables, hydrodynamics variables, or grav¬ 
itational mass as the radial zoning is increased from 
15 to 200 zones. Even for only 10 radial zones the 
error in mass only rises to ~1%. Hence, the zoning in 
the 3D calculations discussed here and in WMM96 is 
probably not a significant source of uncertainty. 

We wish to explore the reliability of the ID model 
as the Lorentz-like factor (Wq — 1) is increased. Fig¬ 
ure H shows the proper central baryon density as a 
function of {Wq — 1). Results are given for two dif¬ 
ferent EOS’s; one for which Mq = 1.64 Mq; and one 
with Mq = 1.70 Mq. Also shown for comparison are 
central densities as a function of average values for 
(IT^ — 1) from 3D calculations. We see that the basic 
trend of increasing central density with increasing or¬ 
bital motion is reproduced, although the central den¬ 
sity in the ID calculations is about 2% higher than 
the 3D calculations for the same average (Wg — 1) fac¬ 
tor. This suggests that replacing the distribution in 
W with a mass weighted average value slightly overes¬ 
timates the effect. Nevertheless, this approach is suf¬ 
ficiently accurate to apply to the schematic parameter 
study of interest here. 

7. Heating 

In WMM96 the released binding energy was as¬ 
sumed to be deposited only in increased fermi energy. 
Any thermal excitation was assumed to be radiated 
away so that the stars remained cold. However, it is 
not necessarily true that the input energy above the 
increasing fermi energy goes into thermal energy or 
that it is efficiently radiated away. If this energy were 
not dissipated, the stars could simply oscillate about 
the equilibrium rather than collapse. We argue, how¬ 
ever, that it seems most likely that such oscillations 


would be quickly damped relative to the time scale for 
inspiraling. Initially, the radial changes will be quite 
small, and the coupling of radial motion to thermal 
excitation could occur, for example, via star quakes 
in analogy with observed pulsar glitches. As the rate 
of energy release becomes more rapid and the crust 
melts, we speculate that the coupling of radial modes 
with the orbital motion, nonradial fluid motion, and 
tidal forces will lead to a complex excitation of higher 
modes and shocks which could further heat the star 
and increase the entropy. Also, the coupling of radial 
modes with the magnetic fields could damp the oscil¬ 
lations. Eventually, as the stars become hot enough, 
T ~ 1 MeV, neutrino viscosity will serve to damp the 
radial motion, but this would be late in the evolution. 

As these dissipative processes come into play it 
seems plausible that significant thermal energy could 
be excited. If the thermal energy is efficiently radiated 
away, then the stars will remain near zero tempera¬ 
ture and the previous calculations are valid. However, 
it is also possible that the energy may not be radiated 
away as rapidly as it is released. In which case the 
damping will be converted into both increased fermi 
energy and thermal energy. An upper limit to the 
temperature of the star would be that correspond¬ 
ing to no radiation during the compression. In the 
present work we estimate the possible heating and 
radiation of the stars as they adjust to the changing 
orbit {Wq — 1) factor and tidal forces. 

If there is sufficient heating, the stars may produce 
an associated neutrino and/or electromagnetic signal 
as they compress. Hence, it is of interest to estimate 
the possible heating of the stars as released binding 
energy is converted to internal energy. 

One can make a simple estimate (Mathews & Wil¬ 
son 1996) for the heating from the change in gravita¬ 
tional binding energy as the stars compress in the 3D 
calculations. From Table |^, as J changes by 4 x 10^° 
cm^ in the three-dimensional numerical calculations 
the angular momentum loss rate is J ^ 1 cm. The 
time to radiate this change in J is AJ/J ~ 1.33 sec. 
From the change of binding energy with central den¬ 
sity for single stars for a particular EOS, it is pos¬ 
sible (Mathews & Wilson 1996) to estimate the en¬ 
ergy available for heating of the stars after increasing 
the fermi energy. This has an associated change in 
baryonic mass, hence it is only an order of magnitude 
estimate which we improve upon below. 

Nevertheless, this change in binding energy could 
correspond to a release of as much as 6 x 10®^ ergs in 


thermal energy and a heating rate of 5 x 10®^ ergs s“^ 
per star. The corresponding average energy over the 
stars could be 2 x 10^® ergs g“^. If this energy were 
injected into a star on the verge of collapse, (central 
density of 2.42 x 10^® g cm“^) it would heat the core 
to a temperature T « 45 MeV (assuming that the 
core has the heat capacity of a degenerate Fermi gas 
of neutrons.) 

This much heating could lead to copious neutrino 
emission and may provide a framework in which to 
produce gamma-ray or X-ray bursts. Hence, there is 
motivation to numerically study this possible heating. 
We do this in the spherical calculations by impos¬ 
ing an energy conserving damping term in the equa¬ 
tions of motion. This damping relaxes the stars to 
their new equilibrium. The damped kinetic energy is 
added as internal energy. By integrating the rate at 
which kinetic energy is damped into Fermi and ther¬ 
mal energy, and estimating the fraction which can be 
subsequently radiated away, we get a measure of the 
possible heating of the star before collapse. 

8. Results 

8.1. Analysis of the Collapse Instability 

With the above ID approximation to the effects 
of orbital motion, we can make a systematic study 
of stellar stability as a function of mass, EOS, and 
(Wg — 1) factor. First we set {W^ — 1) = 0 and run 
a hydrodynamic calculation at zero temperature with 
velocity damping until an equilibrium configuration 
is achieved for a given gravitational mass and EOS. 
Then we increase (Wq — 1) in small increments and 
evolve the star hydrodynamically with conservative 
damping. That is, the damped kinetic energy is added 
to the internal and thermal energy as the calculation 
proceeds. 

For stable stars, we generally observe that the ki¬ 
netic energy at first rises to a maximum and then 
damps to zero in the hydrodynamic simulations. For 
a star which has reached the collapse instability, how¬ 
ever, the kinetic energy first rises to a maximum then 
falls to a minimum and then begins to rise again as the 
collapse ensues. Hence, we define the collapse insta¬ 
bility for this systematic study as the point at which 
the radial kinetic energy begins to increase with time 
rather than relaxing to zero. We wish to analyze the 
heating and neutrino emission up to this point. Once 
the instability is reached, the stars quickly collapse 
and much of the subsequent heating or neutrino emis- 
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sion becomes lost below the event horizon. 

Figures |[^ show the central value of the lapse a 
and the released gravitational energy (in units of 10®^ 
ergs) E 53 as (W^ — 1) increases from zero to the point 
at which dynamical collapse is evident. 

Instability in these calculations is also manifest by 
a decreasing rapidly once it falls below some critical 
value. In these calculations this instability occurs as 
a —> 0.4 — 0.5 depending upon the value of (IF^ — 1), 
the initial mass, and the EOS. A similar phenomenon 
has been noted in previous studies of spherical stars in 
the isotropic gauge. It was noted (Wilson 1979) that 
once a < 0.5 unstable collapse generally ensues, but 
of course those calculations had no — effect. We 

see the same a < 0.5 limit as the mass is increased 
with (W^ — 1) = 0. The central a in the 3D orbit 
calculations near collapse is also about 0.4. 

The size of the Lorentz factor at instability and the 
amount of input thermal energy increases with initial 
gravitational mass Mg and critical mass Mg of the 
EOS as one would expect. The released energy rises 
approximately quadratically with (IFq — I). The total 
energy released up to the instability point increases 
with decreasing mass and increasing Mg of the EOS. 

8.2. Temperature and Neutrino Luminosity 

The absolute surface neutrino luminosity will de¬ 
pend upon details of the neutrino transport from 
the interior. It should scale, however, with surface 
temperature according to a Steffan Boltzmann law, 
Ls oc Tf, and the total luminosity should scale with 
the surface radius r and the lapse Og at the surface, 
Ltot oc L^r'^al . 

If the luminosity is not as great as the rate at which 
thermal energy is added to the stars, then central 
temperatures and the associated neutrino luminosity 
could be quite high. To estimate the luminosity and 
temperature as a function of time consider the New¬ 
tonian angular frequency. 


value ujQ can be estimated by integrating Eq. 

Wo 


w = 


(1 — (64/5)m®/^WQ^^t)^/® 

From Eq. ( p^ ) the time dependence of (m/r) is 
(mwo)^/^ 


m 


^ (I — (64/5)TO^/^WQ'^^t)^/^ 


( 29 ) 


(30) 


Since (IF^ — I) can be thought of as a kind of spe¬ 
cific kinetic energy, it should scale as m/r in a stable 
Keplerian orbit. 


IF^ - 1 = 


YM. 

r 



(31) 


Erom this we get an approximate time history 


(W^o - 1) 


Win 


(1 — (64/5)m 


- 1 


(32) 


where Eq. (^^ is normalized to give the average 
(W^D ~ 1) factor from the 3D calculation which has 
an angular velocity wq. 

Erom this relation of (IFq — 1) as a function of 
time it is possible to now construct a possible pic¬ 
ture of the luminosity as a function of time. Ideally, 
one would like to model the detailed thermal neu¬ 
trino production and transport as released gravita¬ 
tional energy is deposited in the interior. Although 
we have begun such a calculation, a detailed modeling 
of the neutrino transport is quite challenging and it 
will take some time before a systematic study can be 
completed. Nevertheless, we can gain qualitative in¬ 
sight into the signal expected from a simple schematic 
model as follows: 

Erom (Eigs. |[^) we note that the total input ther¬ 
mal energy grows quadratically with (IFq — I), 


A„cx(IF2-I)2 . (33) 


(27) 


and deceleration due to quadrupole radiation (cf. 
Blanchet et al. 1995), 


^ = ^^ 5 / 3 ^ 11/3 
20 


(28) 


The evolution of w from earlier times (t < 0) up to a 


To convert this thermal energy into a luminosity we 
assume that the neutrino flux through any surface at 
radius r should scale as. 


Ey cx 


r^T§ 

PKqT'^ 


d{T^) 

dr ^ p dr ' 


(34) 


where kq is the opacity evaluated at Tq. 

The net neutrino flux passing through r is due to 
neutrinos diffusing throughout the volume interior to 
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r. To approximate the effective temperature of the 
flux passing through r we assume that the thermal 
energy is deposited uniformly in mass throughout the 
star. The neutrino flux at any radius r is then taken 
to be proportional to the mass interior to r, 

Fi, oc m(r) , (35) 


where, 


m{r) =4:TT pr'^dr' . 

Jo 


(36) 


Equations (^) and (^ define an effective neutrino 
temperature profile with radius, (cf. Fig||), 


T(r) = A 


f^m(r')dr'\ /m(r')dr' 


1/2 


(37) 


where A is a normalization constant. We then solve 
equations (|^ through (^ under the boundary con¬ 
dition that T = 0 outside the star, and by equating 
the total thermal energy at any given time to the in¬ 
tegrated thermal internal energy per gram e(T) for 
a given temperature and density profile in our the 
equation of state, i.e. 


e(r)47rr^pdr = Efh 


(38) 


The radial temperature profile of a 1.40 Mq star 
just as the collapse instability is reached, is shown in 
Figure || for the equations of state with Me = 1.70. 
For this EOS central temperatures as high as T w 70 
MeV are possible. 

The neutrino flux is then calculated using Eq. (Q) 
with the proper coefficients included. The flux near 
the surface then gives the luminosity. The luminos¬ 
ity can be used to then define an effective neutrino 
luminosity temperature. 




= 4ttR^ 


4 



(39) 


As the stars compress, released gravitational en¬ 
ergy can be deposited as internal energy to be ra¬ 
diated away by neutrinos. The rate of accumulated 
thermal energy is then given by the balance between 
the rate at which gravitational energy is deposited as 
thermal energy from the contraction of the stars Ein 
and the rate of energy loss L^, by neutrinos. 

We find that the temperature profile as determined 
above is consistent with a scaling Eth ~ where T 


is the evaluated from (Eq. near the the surface. 
This scaling arises because the system can be approx¬ 
imated as a degenerate nucleon gas. Also, from the 
flux scaling (Eq. Q) we find that dT/dr oc T, so that 
L^, ~ T^. Hence, both Eth and L^, scale with the 
surface temperature and we can write 

L.ozT^^kEth ■ (40) 


The evolution of the thermal energy can then be writ¬ 
ten, 

Eth = Ein — kEth ■ (41) 

The constant k is evaluated from Eq. ( ^ above, and 
the rate of deposited thermal energy is evaluated from 
Eqs. ( 1^ and (|3^). 


dE,n d{W^ - 1 ) 

dl/W^ — 1) dt 


(42) 


The analytic solution for the heating of the star then 
becomes. 


Eth 


e 


— kt 


E.ne’^^'dt' 


(43) 


Figure ^ shows the estimated neutrino luminosity 
L„ (in units of 10®^ ergs sec“^), the total accumulated 
internal energy, Eth (in units of 10®^ ergs), and the 
rate of gravitational energy release, Em (in units of 
10®^ ergs sec“^), for a 1.40 Mq star with the Me = 
1.70 EOS. These quantities are plotted as a function 
of time for the last 1.5 sec of a star with Mg = 1.40 
Mq and an EOS with Mp = 1-70 Mq where the time 
scale is defined by Eq. |^. 

The total luminosity can become quite significant 
as the collapse instability is approached. About 4 
sec before collapse, the neutrino luminosity from each 
star rises above 10®^ ergs sec“^. The luminosity ex¬ 
ceeds 10^^ ergs sec“^ about 0.5 sec before collapse. 
The combined neutrino luminosity from the two stars 
ultimately reaches nearly 10®^ erg sec“^ before col¬ 
lapse. This is comparable to the neutrino emission 
from type 11 supernovae, but in this case the emission 
is from bare neutron stars. 

In the quadrupole approximation (cf. Thorne 1980; 
WMM96), the gravity wave luminosity Egw scales 
with the square of the {I + l)th time derivative of the 
mass quadrupole moment, Egw ^ ~ {m/r)^. 

Thus, we write. 


Egw = 




(1 — (64/5)TO®/3tt’g^^f)®/4 


(44) 
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From this we estimate that the power (and angular 
momentum) lost in neutrinos will exceed the energy 
loss in gravity waves for roughly 3 hours before col¬ 
lapse. This means that the late evolution up to col¬ 
lapse may not be determined by the rate of gravita¬ 
tional radiation but by the hydrodynamics and heat 
transport of the compressing neutron stars. 

If the radiative momentum loss dominates at early 
times one may wonder whether it could be observed 
as an increased slow down rate in the orbit of a known 
binary pulsar. We have estimated how much the orbit 
period of the binary pulsar PSR-1913-1-16 would be 
affected. Damour & Taylor (1991) have determined 
that the ratio of the observed orbit period change 
to the general relativistic prediction is p°bs jpGR _ 
1.0081 ±0.0022(galactic) ±0.0076(observational). We 
estimate that any orbit period change from compres- 
sional heating is at least two orders of magnitude be¬ 
low the observational error. 

9. Conclusion 

We have made a survey of the compression, heat¬ 
ing, and collapse of neutron stars in close binaries. In 
particular, we have developed a schematic model to 
describe when the collapse instability may occur as 
a function of initial neutron star mass and the EOS. 
We have also analyzed the possible heating of the neu¬ 
tron star interiors as the stars approach the collapse 
instability. We find that the stars may obtain quite 
high thermal energy and neutrino luminosity in the 
final seconds before collapse. This could have signifi¬ 
cant implications both for gravity wave and neutrino 
astronomy as follows: 

9.1. Implication for Gravity Wave Detectors 

The analysis here indicates that the radiative neu¬ 
trino luminosity could exceed the gravity luminosity 
for hours prior to the collapse instability. If so, this 
could have a profound influence on the inferred grav¬ 
ity wave signal. The loss of orbital angular momen¬ 
tum due to neutrinos and electromagnetic radiation 
will be considerably greater than that of two cold sta¬ 
ble neutron stars. The merger will occur on a shorter 
time scale and the gravity wave signal will be domi¬ 
nated by the dynamics of heating and thermal radi¬ 
ation and not the gravity wave amplitude up to the 
point of the instability. 

Once the collapse instability is reached, we esti¬ 
mate that the formation of one or two black holes will 


occur rather abruptly. After collapse, however, the 
system may not appear simply as two black holes in 
vacuum. As has been observed in supernova calcula¬ 
tions for some time (cf. Mayle & Wilson 1988; Wilson 
& Mayle 1993) this much neutrino radiation is likely 
to ablate electron-positron pairs together with bary- 
onic material from the surface of the stars. Baryons 
ejected in this wind are likely to be present after col¬ 
lapse and may interact with the orbiting objects. To 
some extent they will provide material to accrete onto 
the remaining members (neutron star or black holes) 
of the binary. They may also provide a damping 
medium which could accelerate the decay of the or¬ 
bit. Thirdly, this hot wind material may provide a 
medium in which to anchor the magnetic field lines 
of the precollapse stars (cf. Wilson 1975; RufRni & 
Wilson 1975; Damour et al. 1978). 

We speculate that these effects may serve to accel¬ 
erate the merger of the two black holes. The inter¬ 
action of the stars with this medium may affect the 
dynamics of the black hole inspiral unless the mate¬ 
rial is ejected with sufficiently high velocity. Clearly, 
this is an area which warrants further investigation. 
If our speculation is correct then the gravity wave 
signal becomes a probe of the EOS, hydrodynamics, 
and thermodynamics of the neutrons stars as they 
approach and pass through this collapse instability. 

9.2. Implications for Gamma-Ray Bursts 

The possibility that gamma-ray bursts could be 
generated by neutrino emission from coalescing neu¬ 
tron stars has been speculated upon for some time. 
Recently, Janka & Ruffert (1996) have made post- 
Newtonian hydrodynamics calculations of neutron 
star mergers and included the neutrino emission there¬ 
from. They find high luminosities, but the time scales 
are so short (~ msec) that they conclude that it will 
be difficult to model gamma-ray bursts by neutron 
star mergers. This short time scale stems from the 
time scale for mergers. This difficulty is avoided, how¬ 
ever, in our model in which the time scale is set by the 
gradual compression of the stars. We estimate similar 
luminosities, but in our model the neutrino luminos¬ 
ity endures for much longer times, thus rendering the 
possibility of a gamma-ray burst more plausible. 

We have shown that significant heating and associ¬ 
ated neutrino luminosity is possible in the last seconds 
before the collapse instability. This poses some inter¬ 
esting possibilities for cosmological models of gamma 
ray bursts. The thermal emission of neutrinos pro- 
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vides an environment for the generation of an — e~ 
pair plasma by vi> annihilation around the stars. The 
neutrino emission is occurring in the deepening grav¬ 
itational well of the two stars. Their interactions will 
be enhanced by the curved space around the neutron 
stars. Furthermore, the region between the stars may 
provide an environment for the build up of neutrino 
and matter flux and the production of a pair plasma 
as desired in some gamma-ray burst scenarios (e.g. 
Piran & Shemi 1993). 

In addition to the collapse-induced neutrino emis¬ 
sion itself, the escaping neutrinos are likely to gen¬ 
erate a neutrino heated baryon wind from the stars 
(Mayle & Wilson 1988; Wilson & Mayle 1993). Un¬ 
like in supernovae, the velocity of this wind can be 
be quite high, particularly later in the evolution as 
the neutrino luminosity grows. This later emission 
of the high velocity wind could interact with matter 
emitted previously producing shock heating in envi¬ 
ronments of relatively low optical depth far from the 
stars. The interactions themselves may contribute to 
the production of a pair plasma. 

As a preliminary test of this scenario we ran a cal¬ 
culation of neutron stars instantly heated such that 
the surface temperature was ~ 5 MeV. We then fol¬ 
lowed the neutrino and matter transport using the nu¬ 
merical supernova model of Mayle & Wilson (1993). 
We observed a blow off of the outer layer 10“® 
Mq) of the neutron star. This material was acceler¬ 
ated to a speed corresponding to a relativistic y-factor 
of ~ 10. One possibility is that this high speed mat¬ 
ter interacting with magnetic fields and/or interstellar 
clouds might produce y-rays. 

Finally, we also note that after collapse, the pre¬ 
viously ejected material will continue to experience 
heating either by accretion onto the black holes or by 
ram pressure from the orbiting stars. Once present, 
this plasma might become anchored to magnetic field 
lines around the precollapse stars (Ruflini & Wilson 
1975; Damour et al. 1978). The interactions and 
magnetic recombination of these field lines could also 
contribute to heating and pair plasma production. 

All of these processes may be occurring in the back¬ 
ground of the remaining orbiting binary system from 
times prior to collapse until the final merger to a sin¬ 
gle black hole. This orbit period may lead to an un¬ 
derlying millisecond substructure in associated burst 
signals possibly consistent with observations. 

Clearly, this is an area which also warrants more 


investigation. Work along this line is underway to 
explore such effects as a possible framework in which 
to model cosmological gamma-ray bursts. 
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Fig. 1.— Contours oi - 1 in the Z = 0, X - 
Y plane from the three-dimensional calculations for 
Mg = 1.40 neutron stars with J = 2.6 x 10^^ cm^. 
Contour decrease from a maximum of — 1 = 0.04 
in steps of 0.005. 


Fig. 2.— Central proper baryon density pc for ID 
models (solid lines) as a function of (Wq — 1 ) for a star 
with Mg = 1.40 Mq. Lines are drawn for equations 
of state which give a maximum neutron star mass of 
Mg = 1.70 Mq and 1.64 Mq as labeled. For compar¬ 
ison, the points show numerical results from the 3D 
calculations at various values of (IUq — 1 ) for the Mg 
= 1.70 Mq EOS (squares) and the Mg = 1.70 M© 
EOS (circles). 


Fig. 3.— Central values of the lapse function a and 
the released gravitational energy in units of 10®^ ergs 
E 53 as (Wq — 1 ) increases from zero to the collapse in¬ 
stability point. These calculations use an EOS which 
gives a maximum neutron star mass of Mg = 1-55 
Mq. The different curves are labeled by their associ¬ 
ated gravitational mass. 


Fig. 4.— Same as Figure |3[ but for an EOS which 
gives a maximum neutron star mass of Mg = 1-64 
Mq. 


Fig. 5.— Same as Figure |3[ but for an EOS which 
gives a maximum neutron star mass of Mg = 1-70 
Mq. 


Fig. 6.— Radial temperature profile just as the in¬ 
stability is reached ((lU^ “ 1) = 0.055) of a 1.40 Mq 
star for an EOS which gives a maximum neutron star 
mass of Mg = 1.70 Mq. 


Fig. 7.— Estimated neutrino luminosity (in units 
of 10®^ ergs sec“^), the total accumulated internal 
energy, Eth (in units of 10®^ ergs), and the rate of 
gravitational energy release, Ein (in units of 10®^ ergs 
sec“^), for a 1.40 Mq star with the Me = 1.70 EOS. 
These quantities are plotted as a function of time for 
the last 1.5 sec of a star with Mg = 1.40 Mq and an 
EOS with Mg = 1.70 Mq. 
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Table 1 

Parameters characterizing the orbit calculations for stars with « 1.45 M© and for an EOS 

WITH Me = 1.70. 


J (10^^ CTO^) 

2.2 

2.3 

2.5 

2.6 

2.7 

Mb (Mq) 

1.598 

1.598 

1.620 

1.620 

1.598 

Mg (Mq) 

1.416 

1.420 

1.322 

1.317 

1.423 

/ (Hz) 

410 

310 

280 

250 

267 

dp (km) 

39.4 

40.6 

50.3 

50.7 

53.0 

Pinax (lO^*" g cm“3) 

2.03 

2.70 

3.58 

2.92 

1.93 

1—1 

1 

- 

- 

0.050 

0.052 

0.043 

^min 

0.440 

0.379 

0.283 

0.288 

0.463 

dP 

^max 

1.90 

2.05 

1.78 

2.68 

1.84 

E (M© sec“^) 

0.016 

0.0040 

0.00048 

0.00047 

0.0061 

Orbit 

Unstable 

Stable 

Stable 

Stable 

Stable 

Stars 

Unstable 

Unstable 

Unstable 

Unstable 

Stable 


Note.— Mg is the total gravitational mass of the binary divided by 2. Also, 
/ is the gravity wave frequency, i.e. twice the orbit frequency. 


Table 2 

Parameters characterizing the orbit galculations at the final edit for stars with = 1.40 
(Mb = 1.548) Mq in isolation and for equations of state with Me = 1-70, 1.64, and 1.55 M©. 


J (10^^ cm?) 

EOS - Me 

2.5 

1.70 

2.6 

1.70 

2.5 

1.64 

2.6 

1.64 

2.5 

1.55 

2.6 

1.55 

2.7 

1.55 

Mg (Mq) 

1.300 

1.306 

1.299 

1.306 

1.285 

1.299 

1.310 

/ (Hz) 

270 

287 

272 

251 

256 

260 

233 

dp (km) 

62.1 

64.8 

61.7 

65.6 

55.3 

60.3 

68.2 

Pmax (10^^ g cm“3) 

1.82 

1.79 

1.82 

1.86 

2.95 

2.36 

1.98 

(H^Id - 1) 

0.0384 

0.0365 

0.0387 

0.0354 

0.0434 

0.0298 

0.0339 

^min 

0.435 

0.443 

0.426 

0.444 

0.266 

0.396 

0.440 

4 ^raax 

2.05 

2.02 

2.08 

2.02 

2.91 

2.18 

2.02 

E (Mq sec“^) 

0.0016 

0.00155 

0.00156 

0.00141 

0.000548 

0.00113 

0.00112 

Orbit 

Stable 

Stable 

Stable 

Stable 

Stable 

Stable 

Stable 

Stars 

Stable 

Stable 

Stable 

Stable 

Unstable 

Unstable 

Stable 
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Table 3 

Parameters characterizing the orbit calculations at the final edit for stars with = 1.35 
(Mb = 1.49) M© IN isolation and for equations of state with Me = 1.70 M©. 


J (10^^ cm?) 

1.90 

1.95 

2.00 

2.05 

2.10 

2.15 

Mg (M©) 

1.225 

1.231 

1.238 

1.243 

1.248 

1.252 

/ (Hz) 

413 

388 

365 

346 

330 

312 

dp (km) 

46.4 

48.4 

50.8 

52.8 

54.8 

57.6 

{Win - 1) 

0.0531 

0.0497 

0.0459 

0.0431 

0.0408 

0.0388 

^min 

0.399 

0.413 

0.428 

0.440 

0.449 

0.457 

^max 

2.22 

2.16 

2.10 

2.06 

2.02 

1.99 

E (M© sec“^) 

0.0038 

0.0034 

0.0032 

0.0030 

0.0029 

0.0026 

Orbit 

Unstable 

Stable 

Stable 

Stable 

Stable 

Stable 

Stars 

Unstable 

Stable 

Stable 

Stable 

Stable 

Stable 


Table 4 

Parameters characterizing the orbit calculations at the final edit for stars with = 1.35 
(Mb = 1.49) M© IN isolation and for equations of state with Me = 1.55 M©. 


J (10^^ cm^) 

2.10 

2.15 

2.20 

2.30 

2.40 

2.50 

2.60 

Mg (M©) 

1.238 

1.246 

1.255 

1.262 

1.268 

1.274 

1.279 

/(Hz) 

270 

287 

272 

251 

256 

260 

233 

dp (km) 

50.6 

55.0 

57.8 

62.5 

66.8 

70.1 

74.4 

{Win - 1) 

0.0476 

0.0427 

0.0379 

0.0350 

0.0319 

0.0293 

0.0272 

^min 

0.392 

0.424 

0.452 

0.465 

0.480 

0.493 

0.503 

cfp' 

TYnax 

2.23 

2.10 

2.01 

1.96 

1.91 

1.86 

1.83 

E (M© sec“^) 

0.00022 

0.0023 

0.0022 

0.0019 

0.0016 

0.0013 

0.0011 

Orbit 

Stable 

Stable 

Stable 

Stable 

Stable 

Stable 

Stable 

Stars 

Unstable 

Stable 

Stable 

Stable 

Stable 

Stable 

Stable 
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(A4) 


A. Post Newtonian Analysis of the (W^ — 1) 
Effect 

The effect described in this paper steins from a 
numerical solution to the full Einstein field equations 
in the approximation that the three metric remains 
conformally flat. We have interpreted the apparent 
stellar collapse as the result of a general relativistic 
increase in the strength of the gravitational accelera¬ 
tion which causes the stars to compresses and heat as 
the orbit shrinks. At least in part, this increased grav¬ 
ity involves terms which scale with ~ !)■ Since 
(W^ — 1) can be thought of as a measure of the specific 
kinetic energy, we interpret this part of the increased 
relativistic gravity as arising from the increasing mass 
energy associated with the increasing four velocity the 
binary pair as they approach. 

Although, we expect that post-Newtonian relativ¬ 
ity is a poor approximation in the strong fields near 
the neutron stars, it is nevertheless instructive to look 
for this (W^ — 1) effect in post-Newtonian relativ¬ 
ity. One motivation for such an analysis is that the 
post-Newtonian expansion is an independent approx¬ 
imation to the metric. Hence, it can help to dispel 
concern as to whether the stellar collapse is somehow 
an artifact of our metric choice. This exercise can 
also provide some intuitive insight as to the origin 
of this phenomenon. In this appendix, therefore, we 
outline how a simple post Newtonian expansion can 
exhibit an enhancement of the effective gravitational 
potential which involves terms scaling as (W^ ~ !)■ 

In the post-Newtonian approximation (cf. Wein¬ 
berg 1972) one presumes that it is possible to write 
the metric as the sum of the Minkowski tensor plus 
corrections given by an expansion in powers of ~ 
(GM/r) in a conformally flat metric. For example, 
the time-time component of the metric is written, 

5tt = + 5ti(2) + 5tt(4) • • • , (Al) 


of order and 

r“(0) =p(l-be) 

The remaining T^'^(n) terms then derive from sub¬ 
tracting T‘*(0) from the perfect fluid tensor and re¬ 
taining terms at the appropriate order. 

Let us consider only the gtt component of the met¬ 
ric as an indicator of the strength of the gravitational 
field. For example, the first post-Newtonian correc¬ 
tion to the metric, (?tt(2) term, is just 

V"5„(2) = -87rGr“(0) (A5) 

Thus, we have 

gtt{2) = -2$ (A6) 

where $ is just the Newtonian gravitational potential. 
Velocity dependent terms enter at the next order. Fol¬ 
lowing Weinberg (1972) we have 

gtt{4.) =-2^^ - , (A7) 

where ik is a second gravitational potential given by 

=—+4ttG[T*\2) + T^\2 )] . (A8) 

To obtain insight into dt we first write the appropriate 
terms from the perfect fluid energy momentum tensor, 


T“(2) 


p(l -f e) 



1 ) - 2 $ 


(A9) 


where the approximate identification of (IF^ — I) with 
is valid at order The spatial part is just 

T^\2) = {W^-l)p{l + €) + 3P . (AlO) 

Thus, the source for the Laplacian of includes terms 
which scale as {W^ — 1) times the mass-energy density 
plus a smaller contribution from the second order time 
derivative of the Newtonian potential. 


where the numbers in parentheses (n) denote terms 
of order of u". From this metric, the corresponding 
components of the affine connection and Ricci ten¬ 
sor can be similarly expanded. One also decomposes 
the stress energy tensor into the Newtonian rest mass 
density plus corrections. 

T“ = T“(0)-br“(2)-bT“(4)--- (A2) 

=yb(2)-^Tb(4)... ^ (A3) 

where the numbers in parentheses now denote terms 


2p{l + e){W^ - 1 - ^) + 3P 

(All) 

This suggests how the effective gravitational potential 
might be deeper for binary stars (where (IF^ — 1) >0) 
than the static potential of two isolated stars. The to¬ 
tal compression effect, however, derives from the co¬ 
variant derivative of the stress energy tensor as illus¬ 
trated in Eq. |^. One must consider the effective hy¬ 
drostatic equilibrium of each star which may involve 
terms of higher order than first post Newtonian. 
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